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^S| . We study the effects of repulsive interactions on the critical density for the Bose-Einstein tran- 

sition in a homogeneous dilute gas of bosons. First, we point out that the simple mean field 
approximation produces no change in the critical density, or critical temperature, and discuss the 

ry^ • inadequacies of various contradictory results in the literature. Then, both within the frameworks 

of Ursell operators and of Green's functions, we derive self-consistent equations that include corre- 
lations in the system and predict the change of the critical density. We argue that the dominant 
contribution to this change can be obtained within classical field theory and show that the low- 
est order correction introduced by interactions is linear in the scattering length, a, with a positive 
r^ ■ coefficient. Finally, we calculate this coefficient within various approximations, and compare with 

O ' various recent numerical estimates. 

S: 

I. INTRODUCTION 

4— > . 

'^ , A precise description of the role of interparticle correlations on the Bose-Einstein transition is indispensable to 

understanding its physical nature; indeed, correlations are expected to play an essential role in the very existence of 
superfluidity and related properties, vortices, flow mctastability, etc. In general, dilute systems offer the possibility 
of accurate microscopic treatments. The study of the Bose-Einstein transition in very dilute gases could provide 

'^ , experimental tests of the theory. A large portion of the literature on the modification of the transition temperature is 
^ ■ based on a simple transposition of one of the most popular method of condensed matter physics, mean field theory in 
O , various guises, where the correlations are unmodified by the interactions and remain purely statistical (as in an ideal 
gas). Mean field theories, for instance Gross-Pitaevskii, successfully describe a broad variety of interesting phenomena 
observable in experiments, for example the spatial distribution of the gas in a harmonic trap kM', for a recent review 

\^ ■ of numerous successful applications of mean field theories in Bose-Einstein condensation in atomic gases, see B. Our 
purpose in this paper is to go beyond mean field theories and to explore the effects of correlations on the properties 

fvi ■ of the transition, studying in particular how they modify the transition temperature. 

T-H We assume that the interparticle interaction can be described by a positive scattering length a, equivalent to the 

C interaction of hard spheres of diameter a. We shall also consider a dilute gas, i.e., work in the regime where a is much 

^— ^ smaller than the interparticle distance, an^'^ <C 1, where n is the particle density. The critical number density, nj?, 

^». ' of an ideal gas is given by 
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where A is the thermal wavelength 



n^^X^ ^ C(3/2) ~ 2.612, (1) 



(2) 
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C the Riemann zeta function; m the particle mass and ks Boltzmann's constant. Note that since at the transition 
X ^ n~^'^, the diluteness condition is equivalent to a <C A. 
/\^ ' In an interacting gas, the critical value of the degeneracy parameter, tt-cA^, is modified; the first order change in the 

critical temperature is related to that in the degeneracy parameter by 
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Because mean field theories effectively treat physical systems as ideal gases with modified parameters, the dimension- 
less degeneracy parameter keeps exactly the same value as in an ideal gas.^ To calculate its change as a function of 
the interactions it is necessary to go beyond mean field theories and include correlations arising from the interactions. 

One might expect repulsive interactions to increase the degeneracy parameter - equivalently, to decrease the critical 
temperature, T^, at constant density; in general, the presence of hard cores tends to impede the motion of the 
particles necessary for quantum exchange effects, therefore reducing the influence of quantum statistics and the 
critical temperature. For example, the superfluid transition temperature of liquid *He is below that of an ideal gas of 
the same density. Moreover, applying pressure to liquid ^He effectively increases the role of the repulsive core of the 
potential, and decreases the critical temperature.^ 

Studies on the effect of interactions on the transition began in the 1950's with the work of Huang, Yang, and 
Luttinger [pj , who concluded that the phase transition of the interacting Bose gas "more closely resembles an ordinary 
gas-liquid transition than the Bose-Einstein condensation," but they did not make a specific prediction for the change 
in Tc- Shortly thereafter Lee and Yang Q predicted an increase of Tc proportional to a^/^; later, in Ref. 0] they 
corrected this result and concluded that the shift of the critical temperature is linear in a, with no prediction for 
the magnitude or even the sign of the effect. In 1960, Glassgold et al. predicted again a positive temperature shift 
proportional to a}/'^ Q. Later, Huang predicted an increase ~ a'^/^ [^, and recently ||l^, he predicted that Tc increases 
as a^' ^ , using the same virial expansion as that of Ref. |^ . Despite the lack of qualitative agreement among these many 
solutions of the problem, these studies showed that the changes in question were not merely due to excluded volume 
effects (proportional to the cube of the hard core diameter a) but to more interesting quantum effects, proportional 
to a smaller power of a. 

The problem lay dormant for two decades until it was revisited by Toyoda |ll| , who studied the transition in the 
Bogoliubov approximation in the condensed phase. This work predicted a decrease of the critical temperatures at 
constant density proportional to a^". As the sign agreed with the measurements in liquid ^He, the question appeared 
settled. Toyoda's result was reinforced by numerical Path-Integral Quantum Monte Carlo calculations showing that 
the effect of interparticle repulsion was indeed to decrease the critical temperature |l^ , p^ . Nevertheless, at the time of 
these calculations, the issue did not have the same experimental interest as it has now, and it was not fully appreciated 
that these calculations were limited to relatively high densities and did not explore the region of dilute systems. 

With the prospect of experimental realization of Bose-Einstein condensation in dilute gases, Stoof il4,^ql carried out 
many-body and renormalization group analyses concentrating on the dilute regime. Stoof's work contains interesting 
precursors to the present work, e.g., Ref. |1J| predicts a linear positive shift in the critical temperature about twice 
that of our estimate in mm. Reference |15] predicts more structure in the a dependence of the effect, ^ a In a, p7[ , 
qualitatively similar to the a^ In a we describe below pSi . 

A surprise came when the Monte Carlo calculations for hard sphere bosons were extended to lower densities and 
showed, in addition to the depression of Tc at high densities, the existence of a low density regime where the critical 
temperature is indeed increased by the interaction |l9|l . At very low densities the shift of the critical temperature was 
found to be 
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with c ~ 0.34, determined by a numerical extrapolation to the limit a -^ 0. However, a more recent explicit Monte 
Carlo calculation EQ] of the leading correction to the ideal gas behavior predicts a prefactor c ~ 2.3. One source of the 
discrepancy lies in the non-analytic dependence of Tc on a, discussed below, which gives rise to non-linear corrections 
at the densities where the Monte Carlo calculation of Ref. |l^ was performed. 

In the past several years, the problem was attacked by analytic approaches based on self-consistent non-linear 
equations derived both in the Ursell operator formalism p^, and the Green's function formalism [|l3|. One finds 
in both approaches that the effect of repulsive interactions is to decrease the degeneracy parameter, thus increasing 
the critical temperature at constant density. Moreover, Ref. p6| proves the linearity of ATc in a. This was done 
by observing that the dominant contribution to the shift in the critical density can be calculated by restricting the 
propagators to their zero Matsubara frequency sector, thereby reducing the quantum many-body problem to a classical 



^ Mean field theories can lead to a change of the effective mass W , which in turn affects the value of the thermal wavelength 
in (|l[); with this effect included, the critical value of the degeneracy parameter remains the same as for the ideal gas. 

^In liquid '^He, the repulsive cores similarly reduce the effect of quantum statistics, so that the magnetic susceptibility is 
significantly higher than in an ideal Fermi system with the same density. 



field theoretical problem in three spatial dimensions. An alternative proof of the linearity in a based on renormalization 
group arguments is presented in ||2^. While Refs. pi||l^,p^ all agree on the functional form, they do not provide 
definitive quantitative predictions for the prefactor; Ref. ]2l[| provides c ~ 1, and an estimate in Ref. ||l6[ of an exact 
formula for ATc/T° predicts c - 2-3. In the hmit of a large number N of components H, c = 87r/3C(3/2)4/3 ~ 2.33; 
interestingly, this exact result for N -^ oo agrees with the numerical result c — 2.33 ± 0.25 of Ref. pG| for N = 2. 
The reduction of the problem to classical field theory has been exploited in the recent calculations of the transition 
in classical (j)'^ field theory on the lattice extrapolated to the continuum [E3|24|, which give c ~ 1.3. 

The linearity in a is a non-trivial, non-perturbative result. Since the interaction is itself linear in a, one might 
imagine deriving this result in some form of simple pertubation theory. However, the first order term in a, for fixed 
density, vanishes identically, while all higher order terms have infrared divergences. Nonetheless, various authors have 
attempted to skirt the infrared problems. For example, Ref. [E5[ unjustifiably "regularizes" divergences in sums that 
appear at the transition with an analytic continuation of the Riemann zeta function.R Similarly, Ref. pffl uses a 
virial expansion, unjustified at the critical point, as we discuss below. In another approach, Ref. P6[ attempts to 
exploit differences in first order pertubation theory between the canonical and grand-canonical ensembles in finite 
volume; this pertubative approach necessarily fails in the thermodynamic limit, preventing a direct determination of 
the critical temperature. In Ref. [p7| finite-size-scaling is used to reconcile this approach with the grand-canonical 
calculations. Reference P8| calculates Tc with the help of an "optimized linear delta expansion, which avoids infrared 
divergencies and operates for any N; for N ^ 2 the authors find c ~ 3.0, but the validity of the method is difficult to 
assess, and the accuracy of this result may be affected by uncontrolled errors. 

The aim of this paper is to summarize current understanding of the problem of the transition temperature. We 
provide a more detailed account of our earlier analytical calculations, and in addition compare the Green's function 
and Ursell calculations. The paper is organized as follows: In the next section we recall features of the Bose-Einstein 
transition in an ideal gas and show how the addition of a mean repulsive field does not alter the critical value of 
the degeneracy parameter. Then in Sec. Ill we include correlations and obtain, using alternatively the Ursell and 
Green's function formalisms, simple self-consistent equations which reveal the physical origin of the change in the 
critical temperature. In Sec. IV we show that the dominant contribution to the change in the critical temperature can 
be calculated using a classical field approximation, and we show that the resulting change is linear in the scattering 
length. Section V is devoted to numerical calculations of the coefficient c, and to a numerical exploration of the range 
of validity of the linear behavior. We focus throughout on a spatially uniform system; a discussion of the transition 
temperature of a dilute gas in a trap can be found in Refs. |0,p9^,BO[ . For experimental data in the '^He-Vycor system, 
see|l|. 

II. IDEAL GAS, MEAN FIELD AND RELATED CALCULATIONS 

In a homogeneous system, the number density of a non-condensed ideal Bose gas is given by 

= T393/2[Z) (5) 



(27r)3 eP(^l-t^) _ 1 A3 
with (3 — l/ksT, /i is the chemical potential, z ~ exp(/3/i) the fugacity, and (Jl — 1) 

el = fcV2m; (6) 

the Bose (polylogarithmic) function gp{z) is defined by 
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As /i tends to zero from negative values 33/2(2) —>■ ^(3/2) corresponding to the maximum density for a non-condensed 



gas at a given temperature given by Eq. O) 
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^Indeed, the metliod of Ref. |G3] applied to the simplest case of the non-interacting Bose gas implies that as /i goes to zero 
dn/d/j, = A|n(^(l/2)/r, which is finite and negative, in contradiction to the divergence of the compressibility of the ideal gas. 



The simplest way to include repulsive interactions is in mean field. Assuming that all the effects of interactions can 
be described by an s-wave scattering length, one can generalize Eq. dq) by writing: 
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where the shift of the chemical potential A/i is proportional to the number density: 

PAfj, = 2gnl3 = 4aX^n, 



(8) 
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and g = AttTl a/m. Equation (^ is a simple consequence of the Hartree-Fock approximation, using a pseudopotential 
proportional to a, in which the shift of the single particle energies is given by 



^HF = 2gn; 



(10) 



the factor of two comes from exchange. Since T,hf is independent of momentum we have to increase the chemical 
potential by A/i = S/zf to keep the same particle density as the ideal gas. The same results are obtained in §4 of 
Ref. 111. 




FIG. 1. Plot of the density n\ as a function of the shift of the chemical potential /HA/j,, for different values of the chemical 
potential /3/i. For a given value of Pfi, the self-consistent solution is given by the intersection point (circle) of this curve with 
the straight line /3A/i = iaX^n. 

Because the Hartree-Fock self-energy depends on the density, the relation between the chemical potential at the 
transition and the critical density is more complicated than in the non-interacting case, and the equation fi — A/i = 
is non-linear. Its solution is conveniently obtained with the geometrical method of |^1|], illustrated in Fig. |l|. At fixed 
/i and /3, with fSAjj, as a variable, the density is obtained through (g); then a simple construction provides the value 
of A/i corresponding to the transition. Finally, the density as a function of /i varies as shown in Fig. ^ (full line); 
it behaves similarly to that of the ideal gas. However, the transition now occurs at a positive value of the chemical 
potential and the compressibility (l/n^)dn/dfi is finite, in contrast to the ideal gas where it diverges. As mentioned in 
the introduction, the critical density remains exactly the same as for the ideal gas, because at the transition /i = A/i 
and thus critical density is given by the same integral as for the ideal gas. 
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FIG. 2. Variations of the density nX^ as a function of the chemical potential predicted in mean field (full line), for small 
constant value of the interaction parameter a/A. The dotted line corresponds to the ideal gas. The mean field solution curve 
shows no maximum, but a constantly increasing density as a function of the chemical potential. On the other hand, any finite 
expansion in orders of a/X leads to divergencies around fj, — and may introduce spurious extrema (dashed line). 

It is instructive to use this simple mean field model to test the limit of simple approximations in an exactly soluble 
case. Expanding the right side of Eq. 1^ in powers of a/X we find the virial expansion: 



«A^ =. 93/2(2) -4-5i/2(z) X 53/2(2) + 8 (^-j 33/2(2) |2 [gi/2iz)f + 93/2{z)g-i/2{z)j 
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Novi^, if as Ref. pOJ we consider only the first two terms of the expansion we find that the density (the broken lines 
in Fig. 0) develops a maximum for a negative value of /i; since the density must always be an increasing function of 
/i, it then becomes tempting to infer that a phase transition should take place at this point. As seen in the figure, 
this point corresponds to a smaller density than for the ideal gas; this reasoning would then predict an increase of the 
critical temperature, at constant density, proportional to y/a, precisely the result obtained in |lQ]. But one should 
keep in mind that in this simple model the density maximum is just an artefact of the first order virial expansion, 
as illustrated by the absence of any maximum in the full curve of Fig. g; in fact, inclusion of second order terms 
in a of (|ll|) makes the maximum disappear .Q From the original Eq. (|7|), the critical density cannot change. Similar 
arguments were already given in Refs. ^^ and ||2^; sufficiently close to /i = 0, higher order terms diverge faster and, 
eventually, dominate the lower order terms in any viral expansion.^ 

The simple example above illustrates the dangers of truncating an expansion in a, even within the mean field 
approximation. The physical origin of the difficulty is simple: the very essence of Bose-Einstein condensation is 
the appearance of long exchange cycles over the system, which cluster together all particles that they contain [^ ; 
therefore, the phenomenon is not easily captured within any formalism containing a limitation on the size of clusters; 
further discussion of the effect of long exchange clusters on the position of the Bose-Einstein transition can be found 
in ||19| . One must be very careful in truncating pertubative expansions in which nominally higher-order terms turn 
out to be of comparable order; rather it is necessary in general to sum an infinite number of terms. 

In the calculation of Toyoda ^^, the only mean field taken into account is that due to the condensed particles 
below the transition temperature. His approximation is in fact lowest order Bogoliubov theory. While this theory 



*The second order terms given in Ref. |1C[1 differ from those of Eq. (h_l[); nevertheless they do not change our argument. 

^The end of the discussion of §3.3 of Ref. |B3] was given specifically for the case of attractive interactions; then, instead of a 
density maximum, the naive first order virial correction model predicts the disappearance of the transition, which is replaced 
by a simple crossover between two regimes. For repulsive interactions, the sign of the first order correction is opposite, and the 
density maximum occurs as in Ref. |l(|. 



describes correctly the ground state at zero temperature and its elementary excitations, its extension near the critical 
temperature meets several difficulties; in particular it predicts a first order phase transition g^, a point not taken 
into account. In fact, above Tc, Toyoda's calculation of the free energy is just that of an ideal gas, with no shift in 
the critical temperature. 

III. SELF-CONSISTENT EQUATIONS 

In this section, as well as in the rest of this paper, we concentrate on the non-condensed state and approach the 
critical temperature from above. As we have seen in the previous section, mean field effects which produce merely 
a constant shift in the single particle energies around T^ do not affect the value of the critical temperature. A 
modification of Tc thus requires the inclusion of correlations; whose effect of such correlations is to lower the single 
particle occupation at small k. Thus, as the temperature is lowered, the chemical potential reaches the lowest single 
particle energy for a value smaller than in mean field, resulting in a decrease of jic- 

It is instructive at this stage to consider the highly oversimplified model in which correlations push down only the 
level A; = and all other levels are treated within mean field. Then, at the transition, when the fc = level hits the 
chemical potential, the other particles experience a constant energy shift /i > with respect to the level fc = 0, and 
the critical density is given by: 

Since p, is small, one can expand the Bose function 53/2 (e^'''^) : 



53/2 (e-^'^) = C(3/2) - 2v/^. (13) 

In fact, in calculating the change in the critical density Aric — ric — n!^, one can equivalently expand the statistical 
factor in dl3) at small k: 
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and arrive at the result: 



f°°dkm,^/ 1 1\ 2 ^-— 

We shall often use the approximation ( |l4| ) in the following. We note here that it is valid provided only momenta 
fc ^ A~^ contribute significantly to the integral (^5|), which requires y/2mp, <C A~^. Since we expect the first correction 
beyond mean field to be /3/2 ^ {a/Xf, this condition is satisfied if a <C A. 

As anticipated, the correlations that push down the level fc = lead to a decrease of the critical density, and hence 
to an increase of the critical temperature. Furthermore, the magnitude of the effect is not necessarily analytic in 
the small change p, of the chemical potential, and hence in the interaction strength a. In fact, the expected result 



fifi ~ (a/A)^ together with Eq. (15) lead to Aric/nc ~ a/A. To include correlations more generally we consider two 
approaches, that of Ursell operators used already for this problem in |gl|] , and that of finite temperature field theory. 
Before deriving detailed results, let us spend a moment comparing the two approaches. 

Within finite temperature field theory one typically carries out a systematic expansion of the properties of a many- 
body system, e.g., the pressure, in powers of the interaction V, using either unperturbed Green's functions Gq or 
self-consistent ones, G. Ursell operators provide a different approach to calculating thermodynamic properties of 
an interacting many-particle system, and lead naturally to expansions in terms of correlations of higher and higher 
orders. The Ursell operator of rank n, [/„, describes the correlations of a system of n interacting Boltzmann particles. 
For example, the operator C/27 defined by 

U^ ^ g-/3(p?/2m+pi/2m+y(ri-r2)) _ g-/3(p?/2m+p^/2m) ^ /^gN 

acounts for two-body correlations. One expects that matrix elements of C/„ vanish between states in which one 
of the particles is far away from the others, and, in the tradition of cluster expansions, one writes expansions of 
thermodynamic functions in powers of the Ursell operators C/„. Every term of such an expansions is expected to be 



finite, even for highly singular potentials such as hard spheres. Inclusion of the specific bosonic or fermionic statistics 
gives rise to exchange cycles. 

Solved exactly, both formalisms give in principle identical results for static thermodynamic properties, and detailed 
comparisons of how specific approximations can be formulated in either approach can be found in [ p4[ . In the following 
we study the effects of correlations by means of a simple self-consistent approximation which can be derived in either 
formalism. This simple self-consistent approximation leads to a nonanalytic change in the spectrum at small k. 

A. Ursell operators 

We briefly summarize the principal results obtained in pl| ] with the Ursell method, reformulated here in a way to 
make a ready comparision with the Green's function approach. More details are given in Appendix. 

Reference |2l[] provides the general diagramatic rules to obtain the reduced one-body density operator in momentum 
space, pk, as a function of the ideal gas Bose distribution, fk, and the Ursell operators f/„ (n > 2). Quite generally, 
above the critical point the single particle density operator has the form of a Bose distribution, but with modified 
single particle energies ,|j 

The particle number density is given in terms of pk by 



(2,)""- "'' 

When looking for leading order corrections one can safely ignore Ursell operators, Un, with n > 3. The resulting 
topological structure of the diagrams of the Ursell pertubation series then becomes equivalent to that of the perturba- 
tive expansion using Green's functions with two-body interactions. Furthermore, by treating the matrix elements of 
U2 as momentum independent, one quantitatively recovers the pertubation theory of the Green's function approach 
with a momentum-independent coupling constant related to the s-wave scattering length a. Finally, as we shall be- 
low, at the critical point the density operator p^ becomes large at small momenta, p^ » 1 for /c — > 0, so that the 
approximation 

l + PkC^Pk-[l3iel+Sfik~^i)r' (19) 

can be used systematically (as in (|29| ) below). These remarks explain why the results that we obtain using Ursell 
operators will eventually be identical to those obtained within the Green's function approach in the particular limit 
of a < A. 



^ Since Sfik is real and plays the role of correcting the ideal gas energy in the Bose distribution, the energies k^ /2m + Sfik may 
be regarded as those of statistical quasiparticles, in the sense of Ref. |36| for the Fermi liquid. Such statistical quasiparticles 
are not equivalent to those obtained from the Green's functions. 
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FIG. 3. (a) First order Ursell diagrams; the double line corresponds to one operator U2- Here we slightly change the 
representation of Ref. [21] and replace the dotted segments there (corresponding to a summation over all exchange cycle 
lengths) by closed curves; in this way, the U2 diagrams become very similar to the usual Green's function diagrams. Despite 
this close graphical similarity, the physical interpretation of the diagrams is different: for instance, exchange cycles do not 
appear at all in Green's func tion diagrams, (b) Simplest examples of t he diagrams included in the iteration, (c) Diagrams 
leading to Eqs. ( |24| ) and (106). (d) Bubble diagrams leading to Eq. (109). 



To obtain the mean field result, Eq. (32) of Ref. [£1|, we iterate the first order diagrams shown in Fig. ga (examples 
of iterated diagrams are shown in Fig. gb); this leads to: 



Pk = /fc(Ai- A/i), 



(20) 



with 



(3Afi ~ iaX^ 






■ Pk = 4aA n. 



(21) 



In this approximation the only effect of the interactions is to produce a momentum independent shift of the single 
particle energies which, as discussed in the previous section, can be absorbed in a shift of the chemical potential: 



fj.' = p.- A/i, 



(22) 



leaving the critical density identical to that of the ideal gas. 

To go beyond mean field, we include in the self-consistent equation for p the corrections displayed in Fig. ^c. These 
are formally of second order in a/A and read (§5 of pi|): 



Pk = fk{p' ~Sp.k), 
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with 
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Note how the integral (p3), in which momentum conservation appears explicitly (g is the momentum transfer in a 
binary collision), introduces a fc-dependence of the energy shift, as opposed to the result of the simple mean field 
approximation. 

We assume that the single particle state with fc = still has the lowest energy, so that the phase transition occurs 
when 



/x' - dfik=o = 0. 



(25) 



The critical density is then given by 
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which, because of the fc-dependence of 6fj,k, does not coincide with the critical density of the ideal gas obtained for 
constant Sfj,. Instead 
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In general, Sfik — '5Mfe=o is an increasing function of fc, so that Aric is negative. 
The variable appearing in Eq. (pq), 
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Pk' Pk'-q [Pk+q ~ Pq] , 



(28) 



can be simplified if we notice that when the critical condition (Eq) is fulfilled the dominant contribution to the integrals 
comes from small momenta for which the statitical factors fk diverge. In fact, if the fk^s are evaluated with the free 
particle spectrum, the integral in (pj) becomes logarithmically divergent in the infrared. To see that, we expand the 
fk's as in (O), 



j/3(eO-p'+5pfe) _ ^ 



P{el + Spk - P') ■ 



Setting 

£fe = el + Spk - Spk=Q, 
we obtain the self-consistent relation valid at small fc, 

'aX^V f d^k' f d^q 1 1 
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A simple power counting argument indicates the integral is logarithmically divergent if we replace the self-consistent 
energies e by the free e^. For the self-consistent spectrum, however, no infrared divergences occur, as we shall see in 
III C. 



B. Green's functions 



In the normal state the single particle Green's function G(fc, z^) is given 

G (fc, 2,y) ^ Zy + /i Ej, ^{k, Zjy), 



(32) 



where fc is the single particle momentum, and z^ = 2T:iv / (3 is a Matsubara frequency, with u — Q, ±1, ±2, . . .. The self- 
energy S(fc, z), which describes the effect of the interactions, can be obtained as a series in powers of the interaction 
strength a by standard diagrammatic techniques 0,B^-p9[. The single particle density matrix is related to G by 



Pk 
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G{k,z,) 



(33) 



The criterion for condensation is that the chemical potential p. reaches the bottom of the single particle excitation 
spectrum, and we again assume, as in §111 A, that the lowest single particle state is that with fc = 0. The transition 
point is then determined by the condition ]40|] : 



G-i(0,0) = or E(0,0) = Ai- 



(34) 



At that point. 



G-\k, z,) - z, - el - [E(fc, z,) - E(0, 0)] , 



(35) 



To first order in the interaction strength, the self-energy is given by the Hartree-Fock approxiniation,n leading to a 
contribution 'Shf (see Eq. (10)) independent of both k and z which can then be eliminated by a redefinition of the 
chemical potential, as discussed above. 
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FIG. 4. (a) Green's function diagrams leading to Eq. 
to Eq. (|^, similar to Fig. 3d. 



similar to Fig. 3c. (b) Green's function bubble diagrams leading 



The structure of I](fc,z) in next order is described by the two diagrams in Fig. ^. The second is the exchange 
term of the first, and within the present approximation in which the matrix elements of the interaction do not depend 
on momenta, the two contributions are equal. Replacing the free propagators by their Hartree-Fock version, we have 



E(fc, z^) = 2.g2 



d^k' <fiq f k' jl + f k' +gKl + fk-,)- {I + fk'}fk'+Jk-g 



(36) 



Because the condensation condition ( p4D involves only the Matsubara frequency Zi, — 0, we concentrate from now on 
this contribution. Furthermore, as before, we isolate the dominant contribution by expanding the statistical factors, 
so that fk ~ T/e^, and 



s(M)-s(o,o).-2.^r'/ii;^ 



1 



1 



1 



(27r)3 (27r)3 (^ - ^,'){el,_^ f,') \ £«_, - m' 



^0 _ ,,/ 

'q /-' 



(37) 



Replacing the bare energies e^ by the dressed energies e^ -I- S(fc, 0), observing that /i' = S(0, 0) at the transition, and 



setting 



efc = e° + E(fc,0)-S(0,0), 
we recover Eq. (^. Note that the earlier Sfik — SfJ.k=o is simply E(fc, 0) — S(0, 0). 



(38) 



'^ Implicit in this expression is the summation of two-body collisions via the f-matrix, which relates the two-body potential to 
the scattering length a at low energies. The summation is made with diagrams where the intermediate propagators are free, 
corresponding to two particles interacting in the vacuum plT. 
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C. Discussion 



The change of the critical density is intrinsically related to Siik (or equivalcntly l](fc,0)) by Eq. (^61). In terms of 

C/(fc) = 2m [6fik - Sfik^o] = 2m p(fc, 0) - S(0, 0)] , (39) 



the change of the critical density introduced by the interaction is 

d^k ( 1 1 



Ann = 



(27r)3 1 e'^['=^+^('=)l/2™ - 1 efik^/2m_i 



(40) 



At 






(41) 



where in the second line we have used the approximation (14). 

The heart of the calculation of the critical density is then to determine the function U(k), a non-trivial task, since 
the evaluation of this function by naive perturbation expansion fails because of infrared divergences. However, higher 
order iterations lead to an instability in the energy spectrum at small momenta, as in Ref. |21[| . In the limit of an 
infinite number of iterations, the spectrum around fc = hardens: the self-consistent solution of (|24|) leads indeed to 
£fc ^ fc'^/^, as predicted by Patashinskii and Pokrovskii [Q using the following argument. 

For free particles, the integral in Eq. l pT\ ) contains six powers of momentum in both the numerator and denominator, 
and is thus logarithmically divergent. In order to ensure that the self-consistent solution, Sk, converges in the infrared 
limit, Ek must behave (modulo possible logarithmic corrections) as ~ k" with a < 2, so that the free particle energies, 
el ~ fc2, can be neglected at smah k with respect to S(A:, 0) - E(0, 0). With this behavior, S(fc, 0) - E(0, 0) ~ fc^-^", 
so that we find a self-consistent energy spectrum, Sk ~ S(fc,0) — S(0,0) ~ fc", for a = 3/2. 

The modification of the spectrum occurs only for small momenta k ^ kc, where kc is a scale that will be specified 
below. We note here only that since U{k) is of order a^/A"*, one expects kc to be of order a/A^. For momenta 



oo, perturbation theory becomes applicable leading to U{k)/k'^ -^ 0. The typical momenta involved in the 
0) are of order kc- The validity of Eq. (|4l|) for k <$^ X"^ requires kc <C A~^, which is satisfied in the dilute 



kc<Ck- 

integral 

limit. 

We later present numerical self-consistent solutions of Eq. (|3l]). Here we reconsider the simple analytical model 
calculation of Ref. fl^ which provides an estimate for the scale kc, and acts as a reference for the numerical results 
presented later. In this analytical model we construct a self-consistent energy spectrum at the critical point: 



£fc 



h'k' 

2m 



l](/s,0)- 2(0,0) 



within the approximation (37) for the self-energy, which we write as 

E(M)-mo)^-2/r/|^B„)(jl--l), 

where the bubble diagram contributes 



«<'>-/l75 



d^p 1 

£p+q 



To extract the low momentum structure, below the scale fee, we evaluate the most divergent terms of Eq. 
the following ansatz: 



, n'^k^/^ 



Q{kc-k) + Q{k-kc). 

2m 2m 



With this spectrum, Eq. (Hj) becomes 



B{q) 



Am 



irh'X^kr 



In ( — I -h c 



(42) 

(43) 

(44) 
I using 

(45) 
(46) 
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where c« 2 + 2 In 2 — 7r/2 =1.816, and in the Umit k ^ the self-energy is 



Identifying the right side of this equation with kj h k^^^/2m, the self-consistency condition, Eq. ([l5|), in the hmit 
k ^ appHes that 







^2^^ 


1/2 


a 


k- -- 


= :v. 1 






«20.7— 






U5; 


A2 


A2 



(48) 

As expected, the scale of the low momentum structure is a/A^. However, one should note that the large value of the 
numerical factor implies that the range of validity of the calculation is limited to very small values of a/A (so that 
the condition kc -^ X^^ is fulfilled). 

The energy spectrum obtained with this analytical model is only self-consistent for wavevectors k <ti kc. In the 
limit fc 3> fee we assume that Sk goes over to the free particle spectrum e^ — fi k'^ /2m, ignoring here a logarithmic 
correction (see Sec. V). We smoothly interpolate between these limits, writing 



Thus we estimate the critical temperature as 

ATc 4A f°^ U{k) 



37rC(3/2) 7o fc^ + U{k) 



ian^'^. (50) 



While the precise coefhcient is sensitive to the details of the interpolation between the low and high k limits, e.g., 
([l9|), the result remains of order unity in any case. 

The k^'"^ spectrum is only an approximation and is not stable if higher order corrections are included; from the 
general theory of phase transitions, at T^, Sk ~ fc^^'', where r] = £2/54 ~ 0.02 in an e = 4 — £> expansion or 
•q — 8/{3n^N) ~ 0.14 in the large N limit (with N = 2) |^^. This model provides too strong a modification of the 
spectrum at small momenta. 

The model calculation illustrates however the basic mechanism behind the change of the degeneracy parameter, 
the modification of the single particle energy spectrum at small momentum. This may be understood as a result of 
correlations among particles caused by their repulsive interactions: particles minimize their repulsion by avoiding each 
other in space, i.e., by correlating their positions; the physical origin of the effect is therefore a spatial rearrangement 
that affects the atoms with low momentum. By contrast, atoms with high momenta have too much kinetic energy 
to develop significant correlations. The modification of the spectrum translates into a modification of the population 
of the various levels. In particular the low momentum levels at momentum scale ~ kc are less populated than they 
would in a mean field approximation at the same density, and the overall result is a decrease of the critical density. 

The hardening of the spectrum obtained as a solution of the self-consistent equations, which is responsible for 
the decrease of the critical density, also provides a cure for the infrared divergences which occur in the perturbative 
calculation in second order. However, as we shall see in the next section, such divergences appear in all orders in 
perturbation theory, so that we need a more general scheme to approach the problem. 

IV. CLASSICAL FIELD APPROXIMATION 

In this section we extend the discussion of the previous section in a way that is at the same time more general, in 
that it is not restricted to any particular class of diagrams, and less general, in that only the linear corrections to the 
density are investigated. 

A. Breakdown of perturbation theory 

Our main goal is the calculation of the critical density. As an intermediate step, we distinguish in Eq. (p3) for pk 
the contribution of zero and non-zero Matsubara frequencies: 
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Pk^-TG{k,Qi)~TY,G{k,iLo,). (51) 

The density is obtained by integrating over momentum k (see Eq. (|l8|)). The terms with v ^ are regular at small 
momentum since a non-vanishing Matsubara frequency provides an infrared cutoff. They provide corrections to the 
density, that are analytic in the self-energy, and therefore of the same order as E, starting at order a^ (modulo possible 
logarithmic corrections). On the other hand, the integral for r/ = is singular for small fc, and the infrared divergences 
introduce non-analyticity in a. Since, we are interested here in the dominant correction to the critical density, we 
will retain only this term in pk- Note that the resulting expression for the density is ultraviolet divergent, a problem 
bypassed by calculating the change in the critical density. 

As illustrated by the example of the previous section, infrared divergences also occur in the calculation of self- 
energies E(fc,0); we now use simple power counting arguments to analyze these divergences. Let us first consider 
diagrams in which all the internal lines carry zero Matsubara frequencies. It is convenient here to introduce a new 
notation and set 

el-p' = {e + C^)/2m. (52) 

The quantity ^, a the mean field correlation length, is given by 



T^ 



2mC 



= -p' = -ip-^Hp); (53) 



C plays the role of an infrared cutoff in the integrals. Note that C — > oo (/i' -^ 0) when T ^ T^. In the perturbation 
series, we take the intermediate propagators to be neither free, nor fully self-consistent as in the previous section, 
but containing the mean field contributions. All the functions that are integrated in the diagrams then appear as 
products of fractions of the form 

[K' + CT\ (54) 

where K denotes a generic combination of momenta; it is then natural to use the dimensionless products K(^ as new 
integration variables. Consider then a diagram of order a". The lowest order n — 2 has been already explicitly written 
in (|3^), and it is proportional to (a/A)^ In(fcC), where k is the external momentum. For n > 2, every additional order 
brings in one factor a from the vertex, one integration over three-momenta, a factor T, and two Green's functions 
(the internal lines). The contribution of the diagram can thus be written as: 

where i^ is a dimensionless function, which we do not explicitly need here. The main point is that when one approaches 
the critical temperature, the coherence length becomes large so that the summation of terms ( p5| ) diverges. In the 
critical region, C is '--^ A^/a, so that all the terms in the perturbative expansion are of the same order of magnitude. 
Therefore, at the critical point, perturbation theory is not valid. 

Let us now assume that in a given diagram some propagators carry non-zero Matsubara frequencies so that one 
momentum integration will be altered. For that integration, the presence of an additional imaginary term 2i'Ki'T 
in the denominators of the propagators ensures that no singularity at A: = can take place. Essentially, in the 
corresponding propagators, C, is replaced by a term proportional to A, so that one factor aC,/\^ in (p3) is now replaced 
by a/A. Compared to the diagram with only vanishing Matsubara frequencies, this diagram is down by a factor a/A, 
and thus negligible in a leading order calculation of E. 

B. Classical field approximation 

The diagrams where all Matsubara frequencies vanish are those of an effective theory for static fields. Ignoring 
the non-zero Matsubara frequencies is indeed equivalent to ignoring the (imaginary) time dependence of the field 
operators. In this approximation the many-body problem reduces to a classical field theory in three space dimensions. 

The energy of a classical field configuration is given by 



n^ (fir 



(^ _,>„.„», ^i,<,.,|.)^ (-) 
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The zero Matsubara component of the density is given by (|(^(r)p). By assumption, the wavenumbers of the classical 
field are limited to k less than an ultraviolet cutoff A ~ A~^. As one approaches the critical region, k ^ kc, all the 
terms in the integrand of (|5^) become of the same order of magnitude: 

y, -fc^, (57) 

Am m /i' 

where {T/fj,')k^ is the contribution to the density of the modes with k ^ kc- From Eq. ( |57| ) we see that kc ^ a/X^. 
For k :^ kc perturbation theory in a makes no sense, and in fact all terms in the perturbative expansion are infrared 
divergent. For kc-€.k -^ A~^, perturbation theory is applicable. Note that, in the critical region, C, ~ l/fcc '^ A^/a. 
By a simple rescaling of the fields Lp -^ y/mTip, one can write the effective action for the classical field theory as 

- n/T = jd\ Q|V0(r)p - m^,'\cp{r)\^ + ^I^MI'^) . (58) 

The rescaled fields (j) have the dimensions of an inverse length. The classical theory contains ultraviolet divergences, 
which spoil simple dimensional arguments for the linear change of T^. 

C. Linear dependence of the density correction 

We now consider a diagrammatic expansion of S in terms of the full zero frequency Green's function, defined by: 

-2mG-'^{k) =fc2-2m/i + 2mE(fc,a,G,A), (59) 

from here on we obmit the explicit index j/ = in S and G. In this self consistent expression the self consistent the 
self-energy S depense on /n only through its dependence on G. Instead of /i, we use the dimensionless parameter a 
defined by 

-2m/Li + 2mS(0) =a — . (60) 

A 

The parameter a controls the distance to the critical point; it vanishes exactly at the transition, as opposed to /i. In 
terms of a the Green's function is now given by: 

~2mG-^{k)^k^+a— + U{k), (61) 

where 

C/(fc) = 2TOE(fc, a, G, A) - 2toI](0, a, G, A). (62) 

Since E depends only on the full Green's function, U{k) depends only on a and not on S(0); moreover, the ultraviolet 
divergence in S(fc) is only logarithmic, and the difference U{k) is independent of the cutoff A in the limit A —f oo. 



If we assume that A -^ cx3, the power counting analysis of §IV A implies that 



t^(fc) = ^5(^,«)- (63) 

Inserting this result into ( ^ ) and making the change of variable x ~ kX"^ /a, one finds 

An, = hm / dx L' ^ 64 

nX^a^oJ x'' + a(x,a) 

showing that the change in the critical density is indeed linear in a. 

This result assumes that the limit a-{x, a — > 0) is well defined. This is the case in the self-consistent schemes that 
we discussed above; they avoid the infrared problem of perturbative calculations, and lead to well defined values of 
Auc- Similarly, in calculations involving resummations of bubbles or ladder diagrams the cutoff is provided by an 
effective screening explicitly generated by the infinite resummations. Large N techniques lead to a similar screening, 
with the advantage of also providing an expansion parameter p3,H . 
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On the other hand, situations where the hmit a{x, a ^ 0) is problematic are encountered in perturbation theory 
where, for reasons discussed above, an infrared cutoff is needed; determination of this cutoff through the condensation 
condition can lead to spurious a dependence (an explicit example is worked out in detail in the next section). 

The linearity of the shift in the critical density does not depend on the ultraviolet cutoff and is thus an universal 
quantity. Nevertheless, the universal behavior implicity assumes that the limit A ^ cx) has been taken and is strictly 
valid only in the limit a ^ 0. If a is not sufficiently small, the classical field approximation ceases to be valid and 
non-linear corrections appear. The classical field approximation requires that all momenta involved in the various 
integrations are small in comparison with A ~ A~^ or, in other words, that the integrands are negligibly small for 
momenta k ~ A^^. Only then, for instance, can we use the approximate form of the statistical factors (p9h. This 
requires in particular that k ^ kc <S^ A, yielding a/A <?C 1. In fact, because, as we shall see, the relation between kc 
and a/X involves a large number, this regime is reached only for very small values, a/X < 10~^ — 10~^. 

V. EXPLICIT CALCULATIONS 

Our goal in this section is to provide specific illustrations of the discussions of the previous sections. We first 
present analytical calculations which shed light on the difficulties encountered when attempting to calculate the shift 
in the critical temperature using perturbation theory. Then we show how partial resummations of the perturbative 
expansion generate screening of long range correlations and allow an explicit calculation of the self-energy, and then 
of the transition temperature. Finally we present results of numerical self-consistent calculations, which we compare 
with the analytical counterparts, and evaluate the limitations of the classical field theory. The accuracy of such 
approximative schemes is difficult to gauge a priori. An alternative is to use lattice calculations to solve the three- 
dimensional classical field theory. Results of such calculations have been presented recently |g^,Q . 

A. Second order perturbation theory 

In order to illustrate the difRculties that one meets in perturbative calculation near Tc, let us return for a moment 
to the second order self-energy diagram, which is the lowest order diagram that introduces correlations and therefore 
corrections to the critical density. The value of this diagram for vanishing Matsubara frequencies is given by Eq. (BY 



S(fc) - E(0) = -2g^T^ I 



Sk' d?q 1 



1 



k+q 



fi') eg - M' 



(65) 



(2^)3 (27r)3 (£0,-^0(4-,-/^') 

where g = Ana/m. 

We note that E(fc) is the convolution of three factors of the form I/(fc^ -I- C^^), with ( defined in Eq. (p3|). Using 
the Fourier transform 

''' ^"'^ ^e-^/C, (66) 



(27r)3 fc2 + C-2 47rr 
we obtain 



-r/C 



2mE(fc) = -128^^ i^^j I r'drjoikr) i^^j , (67) 

where jo{x) = sina;/x. This expression contains, as anticipated, a logarithmic divergence at small distances. Let us 
isolate this divergence by separating the Bessel function ^0(2:) into its value at the origin and a correction term: 

3o{x)=jo{0)+(^^~l^ (68) 

The first term gives a momentum independent contribution, given by E(0). Introducing a cutoff 1/A to control the 
ultraviolet divergence, we obtain 



2mS(0) = 128.2 (^)'Ei(- A) « _128.2 ^^y 



ln(^)-7' 



(69) 
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where 7 = 0.577 ... is Euler's constant, and Ei the exponential integral function. The last approximate equality is valid 
when AC ^ 1, which we assume to be the case. The second term, which is regular and equal to 2to(I](A:)— S(0)) = U{k), 
does not require a cutoff. The result is 



3r/C 



C/(fc) - -1287r2 (^) / {j„{kr)-l)^^dr 



= 1287r^ —r 



3 fcC 

-— arctan-— 
kC 3 



I- 






- 1 



(70) 



This equation implies that U{k) is a monotonically increasing function of fc, ^ k"^ at small A:, and growing logarithmi- 
cally at large k. This logarithmic behavior, obtained in perturbation theory, remains in general the dominant behavior 
of U{k) at large fc, i.e., for C""'^ < A;£ A. 

Our result for U{k) can now be used in Eq. ( pd] ) in order to determine the change in the critical density from 
Eq. (^). Because U{k) > 0, this change is negative. We get: 

A 2 r.. ^(fc) 2 p j^(x) 



7rA2 



where we have set 



fcC 



A:c = 87r^ — , 



C/(fc)EEfc2cr(a;), 



J = Cfcc 



(72) 



Note that for small x, a{x) ~ x^ JTItv^ while at large x, (T(a;) ~ 2 (ln(a;/3) — l)/pi^ . The function a{x)/x'^ oc Aefc/e^ 
gives an indication, independent of the specific values of the parameters, of the range of values of x over which 
the single particle spectrum is significantly modified, and hence of the range of momenta contributing to Aric: this 
function monotonically decreases with x, reaching half its maximum value for a; ~ 7, and about 1/10 of its maximum 
when a; ss 20 — 30. Comparison of this momentum scale with the characteristic momentum scale of higher Matsubara 
frequencies (1/A ~ 1/A) gives a constraint on the values of a for which the calculation is meaningful. In particular, a 
has to be small enough that a/ A ^ j^(J/87r^). 

As noted earlier, the momentum dependence of S(fc) is essential for Aric to be non-vanishing. In this second order 
calculation, the (statistical quasiparticle) spectrum remains quadratic at small fc, and is given by 

efc ^ — ( 1 + TT^T^ ) , for k ^ 0. (73) 



2m 



277r2 



As expected, the spectrum of the interacting system is harder than the free spectrum. In the present approximation, 
it is identical to the spectrum of non-interacting particles with an effective mass m* < m. 

The final result for Auc depends on the infrared cutoff Q, which can be determined by the condensation condition 
li' = I](0): 



1 



= ^fc 



In 



3 



■7 



(74) 



In principle the ultraviolet cutoff A could be eliminated by an appropriate counter term calculable in the full theory. 
Alternatively, one could calculate S(0) from the expression (|3^) involving the complete statistical factors. The result 
of such a calculation would be to replace the term In(AC) in Eq. (^) by ln(C/A), up to a numerical additive constant. 
Here we shall simply choose an ultraviolet cutoff A = 1/A, keeping in mind that there is arbitrariness in the procedure 
which affects the final result, since for a/A ^ 1, then C ^ A, and ln(C/A) will eventually dominate. Nevertheless, it is 
intructive to solve the equation above for ^ as a function of a with this choice of cutoff. Typical values are given in 
the table below: 



a/A 


C/A 


J 


c 


0.01 


6.5 


5.1 


2.4 


0.001 


23 


1.8 


0.90 


10-^ 


153 


1.2 


0.60 


10-^ 


1208 


0.95 


0.47 










10-*^ 


7.510^ 


0.59 


0.29 


10-^" 


7.010^ 


0.55 


0.27 


10-^^ 


6.510* 


0.51 


0.25 
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The behavior of ( with a is understandable: if a is large, condensation takes place far from the mean field value, hence 
the small C- If a. is small, condensation takes place near the mean field value for which ( — > oo. In fact Eg. ([74| ) shows 
that C "^ y? ja up to a logarithmic correction. The last column of the table gives the coefficient c in Eq. ^ for ATc- 
The variation of ATc with a follows closely that oi J J dx a{x)/x'^ = J/Gtt; that is, the term in J^ in the denominator 
plays almost no rolejj 

This simple calculation also illustrates the limits of a pertubative approach. The infrared cutoff introduces a new 
scale in the problem which spoils the argument leading to the linearity of the a-dependence of ATc (c is not a 
constant). The condensation condition ( [7^ ) which relates the infrared cutoff to the microscopic length A, induces a 
spurious logarithmic correction which does not vanish as a — » 0. 

B. Non self-consistent bubble sums 

The previous calculation illustrates how the mixing of ultraviolet and infrared divergences in perturbation theory 
can produce spurious a dependences. It is therefore desirable to find approximations in which the infrared cutoff 
is internally generated. One such approximation was already presented in Sec. IIIC. We turn now to another, the 
resummation of bubble diagrams, as illustrated in Fig. 0b. Again, the quality of such an approximation can only be 
gauged by a comparison with an exact calculation, except in large N limit where the bubble summation becomes 
exact itself ||,||l. 

The one bubble diagram can be calculated explicitly. Keeping an infrared cutoff, we have 

, , f (fip 1 47r gC 

""^'^^^J (^ {el - l/2mC^)(eO,, - 1/2^^) ^ A^ ^'^^^^^T' ^''^ 

In the infinite cutoff limit {C, — > oo) this simplifies into 

The Fourier transform of B{q) is nothing but the leading contribution to the density-density correlation function. At 
the critical point this correlations behaves as 1/r^, so that density fluctuations are correlated over very large distances; 
this is the physical origin of the infrared divergences of perturbative calculations. Nevertheless, these fluctuations can 
be screened, for instance by summing the bubble or ladder diagrams. The respective contributions of the two classes 
of diagrams actually differ only by the number of exchange diagrams. For the bubble sum, the correlation function 
reads 



l + 2gB{q) X^Tq + k^ 
and is now regular at small q. The screening wave number fee is given by 

In the ladder approximation the factor of 2 in the denominator of (O) is absent, and, correspondingly, fee — 47r^a/A^. 
Now, an infrared cutoff is no longer needed in the calculation of U{k), and the limit C, ^ oo can be taken. One 
finds 



2 1'°° I 

U{k) = ^rkln / dx 

TT^ Jo 1 + = 



^In 
2 



1 + a; 



1-a; 



(79) 



where k = k/kc- To study the limiting behavior of U{k) at large and small fc, it is convenient to transform this 
expression as follows. First we integrate twice by parts to obtain 



*Note that these second order results are closely related to the pertubative calculation of [El| where c ~ 1 was obtained by 
looking at values 0.001 < a/A < 0.01. 
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2 fc2 r°° 2{l + xk)[\ii{1 + xk) - 1] 






Taking the derivative of the integrand with respect to k which obtain 

d r°° , 2(l + a;K)nn(l + XK) - 1] 

' ax ■ 



dn Jq (1 — x-^)^ 

= — dx ln(l + xk) I — -7- — — -7- I = In k. (81) 

2 J ^ '\{l-xy (l + xyj 1-k2 ^ ' 

The K integral can now be expressed in terms of the polylogarithmic function gp{x), defined in Eq.(|^); the integration 
constant is choosen to make U{k = 0) = 0. Thus, 



^ ^[1 _ InK] + -hiKln(l + k) + - [52(1 - k) +52(-k)] - -„ 

TT^ K Z 2 12 



U{k) = - — - \ 4^ -1^^] + :^ ln«ln(l + n) + - [52(1 - «) + 52(-aj)] - — h (82) 



However, to derive the hmiting behavior of U{k) it is more convenient to take the hmits in Eq. (|8l| ) and integrate 
afterwards. 

For small k (k = k/kc ^ 0.1), U{k) is well approximated by its small k behavior: 

^(^) = " 5^^' (1^ £ " ' ^ « ^- (^3) 

As expected from perturbation theory, U{k) grows logarithmically for large momenta, fc, and for k/kc <^ 50, U{k) is 
well approximated by: 

2fc? /. k 



TT 



fee 



[/(fc) = ^ In — - 1 , k/kc > 1. (84) 



From the small k behavior of fc^ + U{k) one can estimate the critical index 77. The logarithmic term indicates 
a modified power law in the low momentum limit ^ k^^^' ~ k^{l — rjlnk + ...). Comparing the coefficients of the 
logarithmic terms we obtain 

'^ = ^ ~ 0-068. (85) 

Due to the exchange contributions this value differs by a factor of 2 from the usual large N resultsjj 
The change in the critical density is now 

Let us first estimate the range of kq — fco/^c where a{K) dominates over k'^. Using the small k asymptotics of U{k) we 
estimate kq ^ exp[— 37r2/2] ^ 1. Therefore, we can again ignore the term in a{K) in the denominator of ( pq ) without 
making a significant error; it only brings in an harmless singularity at small n. We get then: 



4:krX f°° dn f°° dx fx 

— In 



7r3C(3/2) Jo K Jo 1 + XK\2 



l + x 



1-x 



(87) 



®Note however that the expansion in powers of rj is meaningful only if the magnitude of i] is controlled by a small parameter, 
such as in the e-expansion or the l/A''-expansion. The estimate presented here should therefore not be viewed as a particular 
prediction for the critical index 77; it gives nevertheless an indication of how the spectrum is modified at small k by the 
resummation of particle-hole bubbles. Another estimate of the effect of bubble summation was presented in Ref. ( [L6[); there 
we tried to estimate the change of the spectrum with respect to the fc^" self-consistent solution. Once the bubble sum is 
included however, self-consistency does not further alter the spectrum at low momentum, as later in this section. As a result, 
the exponent rj that one finds here is much smaller than the crude estimate in Ref. ( |_6|). 
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In order to calculate the integral, we want to exchange the orders of the k and x integrals. Since the integrals, however, 
are not absolutely convergent, before we do so we need to introduce a regularization, inserting a factor k'^ in the n 
integral, and taking the limit e -^ 0+. With this factor we may exchange the orders of integration. The k integral 
becomes 



, k"-^ 1 

an = 

Q 1 + XK ex'^ 



The remaining x integral becomes 



dx ; 



^1 
2^" 



1 



1 



For e = this integral vanishes identically. Thus we may replace x 
remaining integral is 



(89) 



by a; "^ — 1 which goes to — elnx as e -^ 0. The 



da; In a; — In 
2 



1 



l-x 



-1 =-^. 



The factors of e cancel out, and we find 



Ar 



krX 



Ait 



Uc 2^(3/2) C(3/2)A' 

We finally obtain the changes in the transition density and the transition temperature: 

Aric 47r .,. ATc 2 An^ 

= ^nri ' — 

3 ric 



C(3/2)4/3 



Tc 



2.33 anl/^ 



(90) 



(91) 



(92) 



This result for the bubble sum agrees with the leading order result of the 1/N expansion. It is interesting to observe 
that the leading order 1/N result is independent of A^. Since aN is kept constant in the 1/A^ expansion, kc is effectively 
independent of N {kc — 27r^aiV/A^), while U{k) is of order 1/A^. Therefore, the approximation of neclecting U{k) in 
the denominator of Eq. ( p6|) is justified in the 1/N expansion. 

In the bubble sum, we can keep U{k) in the denominator and calculate the integral in Eq. (|8^) numerically, and 
find a reduction the linear coefficient of the critical temperature shift from c = 2.33 to c = 2.20 for N = 2. In this 
approximation the condensation condition reads 



2mE(0) = §■ 



A/ka 



dx- 



2k kc 
^^In- 

7r2 A + 



which gives the mean field correlation length. 



1 

e 



2fc2 



In 



Xkc 



(93) 



(94) 



As before we have taken A — 1/A and assumed that Xkc ^ 1. As opposed to the second order calculation, Eq. (|7J), 
the condition (M) does not mix the infrared and ultraviolet cutoff, and does not introduce any spurious a dependence 
in the final result for the shift in the critical temparature. 

The condensation condition ( |93| ) is the only place where the microscopic scale A enters explicitly. However, the 
classical field theory result for ATc assumes implicitly that the contributions of momenta k ~ A~^ are vanishingly 
small. Alternatively, if we were to cut the integration in ( p6|) off at fc ~ A, one should find a result independent of the 
specific value of A ^ A~^. In fact we have seen that the momenta important in the determination of ATc are k ~ kc. 
The validity of the classical field approximation requires that kc ^ A^^, or, since kc ~ Stt^oA^, a/A ^ 1 
the linear regime is attained only for anomalously small a. When a is not so small, non-linear corrections ~ a 
appear, which tend to decrease the value of ATc, as discussed in Ref. |18[ (see also below). 

In Ref. PI, Stoof examines the appearance of Bose-Einstein condensation, calculating the shift of the critical 
density within a real time formalism. The approach includes, not only the mean field contributions, but also sums of 
ladder graphs within the many-body T-matrix-approximation. He derives an analytical formula for the modification 
of the energy spectrum, from which Stoof obtains a relative increase of the critical temperature, 4.66 an^'^, exactly 
twice the value of the large-N calculation |23]. Summing ladders, and neclecting U{k) in the denominator of Eq. ( pq ) 
we indeed reproduce this result. Evaluating the entire integral numerically, one obtains c = 3.90. 



TT^; thus, 
2 ln(a/A) 
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C. Self-consistent calculations 



We now solve numerically the self-consistent calculations discussed in Sec. IIIC. We quantitatively compare three 
different approximations for the self-energy, the one bubble approximation, Eq. ([43|), 



S(fc) - I](0) 



-2ff'r 



(27r)3 



Biq) 



1 



1 



£k- 



and, to compare with previous calculations, the ladder summation of particle-particle scattering processes 



S(fc) - I](0) = 



-2g^T 



dh 



Biq) 



(2^)3 1 + .gB(<7) \ek 



1 



1 



and, finally, the bubble summation of particle-hole scattering processes 

d^q B{q) ( 1 



E(fc)-E(0) = -2,^T/^- 



- 2ffi?(<z) 



^k — q 



(95) 



(96) 



(97) 



The energy spectrum ej, in the denominators are determined self-consistently using Eq. (H3); B{q) is given in Eq. (H). 
Although the integrals in Eqs. (P5|)-(p6|) giving the difference of the self-energies, U{k) — 2m[E(fc) — S(0)], are 
convergent, we introduce a large momentum cutoff A for their numerical evaluation {U{k > A) = 0). Only in the 
limiting case A — > oo, will U{k) become independent of A; for any finite cutoff, the energy spectrum depends weakly 
on A. The cutoff enters only through the dimensionless parameter AA^/a. For the numerical calculation the value 
AA^/a ~ 800 was used. For the self-consistent bubble calculation we further studied the influence of the cutoff to 
extrapolate numerically to the limit A — > oo. 
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FIG. 5. Self-energy U(k) in units of a^/A"* plotted as a function of a; = kX^ /a for the three approximations discussed in the 
text: (a) self-consistent one bubble approximation, (b) self-consistent ladder summation, and (c) self-consistent bubble sum. 
These result are obtained with an ultraviolet cutoff of AA^/a ~ 800. The dashed line shows the analytical (not self-consistent) 
calculation of U{k) in the bubble approximation, Eq. (p2|). 
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Figure || summarizes the numerical results of this section in terms of the self-energies S(fc) — E(0) corresponding 
to the three different approximations. The various curves in Fig. display the logarithmic growth of U{k) at large k. 
Note however that within the present approximations the overall magnitude of U{k) is determined by the behavior 
of the spectrum at small k: the harder the spectrum, the larger U{k), and the larger the value of c. Nevertheless, the 
values of the shifts in the critical temperature remain comparable. For instance, for the value of the cutoff given above, 
the shifts of the critical temperature that we obtain from Eq. mVt) are: ATc/T^ ~ 3.8 an^/'^ for the self-consistent 
one bubble calculation, ATc/T^ ~ 1.6 an^/^ for the self-consistent bubble sum, and ATc/T° ~ 2.5 an^/^ for the 
self-consistent ladder sum. These values still depend weakly on the value of the ultraviolet cutoff A, and still conatin 
logarithmic corrections ~ a^ ln(a/A), as we shall see below. 

We now compare in more detail these results with our analytical calculations and discuss briefly the extrapolation 
A — > oo, i.e. the extrapolated result of the bubble summation is ATc/Tc ~ 2.0an^/^. 

1. Self- consistent one-bubble calculation 

In the limit fc -^ 0, we expect to recover the k^^"^ behavior of the analytical model of Sec. IIIC. By fitting the 
numerical data to the following functional form 

l,l/2p/2 

2m(S(fc) ~ I](0)) 



1 + ai/2(fc/A;e)i/2 + a,{k/k,) + ay^ik/kc)^^ + ... ' 

we extract a momentum scale kc, which agrees quantitativly with that of the analytical calculation, kc ^ 20 a/ A^. 
However, the spectrum very soon deviates from this behavior, due to the large value of the coefficient ai/2 — 0.9. 
At intermediate wavevectors, around fee, U{k) is roughly linear, and eventually grows logarithmically for k ^ kc, as 
expected from pertubation theory. 

2. Self- consistent bubble sum 

As we have seen in a previous example the main effect of self-consistency is to modify the spectrum at low mo- 
mentum, avoiding infrared divergencies. Since, however, the bubble sum already provides a screening of the long 
range correlations leading to the infrared divergences, we do not expect qualitative changes in U{k) in going from 
the non self-consistent result of Sec. VB, to the fully self-consistent calculations. This behavior can be seen in Fig. |^: 
deviations occur only at high momenta, mainly due to the influence of the finite cutoff in the numerical solution. 

To study more quantitatively the infiucnce of a large but finite cutoff A we have performed a self-consistent calcu- 
lation of U{k) numerically for several values of AA^/a. As explained in Rcf. iQ we expect a logarithmic dependence 
on K\^ / a; therefore we have used this functional form to fit our numerical data, which provides 



— ^~ 1.95an^/3 



1 + 32-r^ln 21— — 
AA V AA 



A/a > 1. (99) 



Extrapolating to A ^ oo we obtain c ~ 2.0, which is slightly smaller than the shift obtained for the non self consistent 
bubble sum. Alternatevely, taking a finite value for AA that is independent of a, e.g. AA ^ 1, provides a logarithmic 
correction which limits the linear regime to very small values of av}'^ . The precise value of this correction is model 
dependent, as we see in the following subsection. 

3. Influence of non-zero Matsubara frequencies 

Non-linear corrections to the critical temperature shift cannot be obtained within the zero Matsubara frequency 
sector. One possibility would be to use an effective field theory which includes the effects of non-zero frequencies, 
e.g., as new vertices in the effective action. Here, we use a different approach, solving the following pair of non-linear 
equations, 

., = ^ + S(fc)-E(0), S(fc)-E(0) = -2.g2r/^ J^nf^--)' (^OO) 

2m J [2ny l + 2gB[q) \ek-q eq J 

21 



where the cutoff in the bubble diagram integral 



B{q)-Pj J^fpfp+q 



(101) 



is no longer a simple step function, but rather the smoother Bose function fk — (exp(/3efe) — 1)"^. We calculate the 
shift in the transition temperature using Eq. (Eoh. Although this does not correspond to a systematic approximation, 
it provides an illustration of the effect of keeping the full statistical factors in the calculation (instead of using their 
classical limit). 
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FIG. 6. Dependence of t he t ransiti on t emperature, Tc/T^, of a dilute homogeneous Bose gas on scattering length a/A 
calculated by solving Eqs. (IOC) and (101) self-consistently. The dashed line is a fit to the data points, given by Eq. (102) 



The asymptotic linear behavior extracted from this fit is shown for comparision. The linear behavior is seen only at very small 
values of a/ A. 

In Fig. 1^ we show the calculated critical temperature in the dilute region. The Bose functi ons in (|4(]| ) lead to a^ In a 
corrections in l^Tc/Tc. On the other hand, the Bose functions in the bubble diagram, Eq. (101) lead to less singular 
corrections. We ignore them here and fit the the numerical datas to the same functional form of Eq. (|99| ) as in the 
last subsection, and find, 



at; 



l.Qan^/^ \\ + 2.&an^/^ \n{'i.lan^'^) 



(102) 



Even in the very dilute region, n^'^a ~ 0.01, the logarithmic corrections are noticeable and reduce the temperature 
shift with respect to the linear prediction. This provides a possible explanation for the discrepancy of the different 
Monte Carlo resul ts I2(]|p^ , p4| and [^; whereas Refs. [ P0|j2^j24[] calculated the linear corrections directly in the limit 
v}/^a -^ 0, Rcf. [|l9[ performed several calculations in the density regime 10~^ < na' 
critical density much smaller than expected from the linear formula of Refs. P0| , P3|,P4| 



^^,0.1 finding a shift of the 
Although the logarithmic 



corrections tend to decrease this linear shift, the approximations underlying Eqs. (100) and (101) are too crude to 
allow quantitative comparision. 

In p5| Bijlsma and Stoof, using renormalization group techniques, obtained an increase of the critical temperature. 
A peculiar feature of their results is that the dependence of the critical temperature on the dimensionless parameter 



an 



1/3 



is given by an unusual curve, going as a In a in the limit of vanishing interaction |17]|. The interpretation of such 



an unexpected dependence is not clear at this stage. 



VI. CONCLUSION 

In this paper, we have studied the effects of particle interactions and correlations on the transition temperature 
for Bose Einstein condensation, and derived the leading effects beyond mean field in dilute systems. Our study is 
general and not limited to any particular approximation, for instance an arbitrary selection of class of diagrams in a 
perturbation expansion. We have shown that the leading term in the change of the critical density is first order in the 
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scattering length a, and can be derived by solving the corresponding classical field theory. Estimating analytically the 
coefficient requires in general uncontrolled approximations. Among the various approximations that we have tried, our 
prefered result is the self-consistent calculation of the sum of bubble diagrams, which gives a coefficient 2.0. However 
this number should not be trusted at a 10% level. Compared with the most recent numerical results of Refs. p3| , p4| , 
c ~ 1.3, our value is still acceptable; the complexity of the mathematical problem does not permit one to make a 
definitive prediction of the prefactor of the linear term from an analytic analysis. 

It is remarkable, that, despite this complexity, all approximations that we have used lead to comparable results: to 
get the right order of magnitude of the critical density or temperature change a precise determination of the energy 
shift U{k) is not required. The contribution of this function to ATc are actually close to "all or nothing" for extreme k 
values: for small fc, the function is larger than the free particle energy and the corresponding momenta are completely 
depopulated, the precise value of U{k) is not relevant; for large fc, the free particle energy dominates and the value of 
U{k) is also irrelevant. The important feature of U{k) are the crossover values at which it is comparable to the free 
particle energy spectrum, and the way this region is crossed by the function. 

We have limited ourselves to an homogeneous gas contained in a box, ignoring the influence of a possible external 
potential, for example magnetic traps and optical lattices. In both such systems, the dimensionality can vary con- 
tinuously from three to two or smaller, and, therefore, affects the nature of the transition. We will discuss them in 
future publications. 
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APPENDIX 



Since the formalism of Ursell operators is less common than that of Green's functions, we give in this appendix a 
few more technical details concerning the equations written in § III A; this will allow the interested reader to make 
contact with the calculations of ref. pl| more easily. For instance, the right side of Eq. ( pi] ) can be obtained from 
Eqs. (58), (53) and (55) of this reference, which provide: 



X,, 



e'^'p,. 



(103) 



so that (55) becomes: 



log 



4a 
T 



d^k pk X e^^^f" 



(104) 



Equation ( |21| ) is then nothing but the first order term in an expansion of this result in powers of a. 

Similarly, Eq. (gj) can be obtained as the lowest order expansion of a relation obtained from Eqs. (55) and (81) of 

ill, 
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log 



(«)%3/5A^j^ 



(105) 



with 



J2 = A« 



d^k' 



d^q 



{2nrJ (2^) 



-3 Pk' Pk'-qPk+q X e 



PA(k,k',q) 



(106) 



and 



A{k,k',q) == 6pk' +Spk 



5p, 



k+q- 



(107) 



We note that we have changed the sign convention of pit] by introducing a minus sign in the right hand side of (EJ); 
in this way, positive Ajj, as well as positive 6pk correspond to positive corrections to the self-energies. This convention 
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makes more straightforward the comparison between 5^k and the self-energy E(fc) introduced in the Green's function 
formahsm. 



The exact form of the 5/^'s is not important for the discussion of § III A. In the context of mean field, what 
matters actually is only the existence of some k independent form of A/x, and one could use expression (104) as well; 
nevertheless, it would not improve the accuracy either, since it , but it is actually just a consequence of the simplest 
approximation used for the self-consistent equation for p. As for correlations effects, the only essential property is the 
momentum conservation rule that appears in (E4h as well as in ( |106| ). 

Higher orders can be readily incorporated into the self-consistent equation pl| ; for instance, a summation of bubble 
diagrams shown in Fig. 0d leads to the generalization of Eq. 



(i{5iXk - Spio) = -4^A^ / ,f ^3 ^ "^y, , (Pk+g - P,) (108) 



A J (27r)3l + A(g) 
where 



A{q)^2jX' J^Pk'Pk'-,. (109) 



With the factor A{q) in the denominator the integral of Eq. ( 108) is convergent in the infrared with a free particle 
spectrum. Further generalizations are discussed in pl|. Numerical solutions of particular approximations are presented 



in Sec. V. As far as the bubble summation of Eq. (108) is concerned, we remark that a summation of Ursell ladder-like 



diagrams leads to the same result without the 2 in the denominator; for more details, see p4 
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